Effects of random localizing events on matter waves: formalism and examples 
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A formalism is introduced to describe a number of physical processes that may break down the 
coherence of a matter wave over a characteristic length scale I. In a second-quantized description, 
an appropriate master equation for a set of bosonic "modes" (such as atoms in a lattice, in a 
tight-binding approximation) is derived. Two kinds of "localizing processes" are discussed in some 
detail and shown to lead to master equations of this general form: spontaneous emission, and 
modulation by external random potentials. Some of the dynamical consequences of these processes 
are considered: in particular, it is shown that they generically lead to a damping of the motion of 
the matter-wave currents, and may also cause a "flattening" of the density distribution of a trapped 
condensate at rest. 



I. INTRODUCTION 

Experiments in which the wavefunction of a material 
particle (or system of particles) becomes delocalized, in 
a coherent way, over a relatively large region of space 
have become increasingly common in recent years, even- 
tually giving rise to the entire subdiscipline of matter- 
wave physics. The purpose of this paper is to charac- 
terize, and study some of the consequences of a num- 
ber of physical process that may break this long-range 
coherence, in particular, for large numbers of particles 
described in a second-quantized formalism. These (as- 
sumed random) events will be referred to here as "local- 
izing events" because, even though they do not necessar- 
ily pinpoint the actual position of the particle (as long as 
they go unobserved), the breakdown of coherence over a 
scale I results in a density operator that can be consis- 
tently interpreted as describing an ensemble of particles 
localized in regions of space of a size ~ /, centered at 
some unknown point (s) xq, with the xq ranging over the 
region occupied by the original coherent wavefunction. 

Traditionally, the study of these types of processes has 
taken place mostly in the context of research into deco- 
herence and the quantum-to-classical transition (see 
for an early review, and 

2] 

for more recent work and 
further references). Interferometric experimental obser- 
vations of decoherence due to particle scattering Q and 
emission of thermal radiation |4| have recently been re- 
ported; these papers also include many references to ear- 
lier relevant experiments, such as, e.g., [H, |(| Q • 

In this paper, a number of physical processes that 
may result in "localizing events" will be described, fo- 
cusing in particular on systems of cold atoms (including 
Bose-Einstein condensates) in optical lattices. A simple 
master-equation formalism will then be introduced to de- 
scribe the decoherence induced by these processes in a 
second-quantized picture, and it will be shown that their 
effect on the transport properties of the system takes 
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the form of a damping of the matter- wave current. In 
one particular case, namely, that of modulation by ex- 
ternal random potentials, an expression for the damping 
is derived that is very similar to one recently proposed 
Q to explain the damping observed in certain one- 
dimensional atomic Bosc-Einstcin condensates. Finally, 
the effects of these random localizing events on the den- 
sity distribution of a trapped condensate at rest are also 
briefly discussed. 

It may perhaps need to be stressed that words such 
as "localization" or "localizing" are used in this paper 
in the sense described above (i.e., breakdown of long- 
range coherence), and not in the technical sense of, e.g., 
the Anderson localization transition, or the so-called "dy- 
namical localization" phenomenon. As regards the first 
of these, the processes considered here are all dynamic 
in nature, whereas Anderson localization |10| . in its sim- 
plest form, concerns the nature of the eigenfunctions of 
a particle in a disordered but static potential. As for dy- 
namical localization, it is a phenomenon that takes place 
specifically in the quantized versions of classically chaotic 
systems, which does not really describe most of the sys- 
tems to be considered here. That being said, however, 
it should also be pointed out that some of the systems 
described in Section V below may, under some circum- 
stances, exhibit classical chaos and/or Anderson localiza- 
tion [TJ , so it might be interesting to study the interrela- 
tionships between all these phenomena for such systems. 
This may be pursued at some later date. 



II. FORMALISM 

A. Localization by direct wavefunction collapse 

Suppose that we have a single particle in a state de- 
scribed by a wavefunction ip(x), which we imagine as 
extending over a relatively large region of space. A "lo- 
calizing event" is one that would cause the wavefunction 
to collapse to a region of space of size I (we work in one 
dimension throughout) around a point Xq; the resulting 
wavepacket may then be described by a wavefunction 
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f(x — Xq), where the function f(x) has width I and is 
centered at x = 0. Clearly, we expect the overall prob- 
ability of such a collapse to occur to be proportional to 



p(xo,f) 



f(x — xo)*ip(x) dx 



(1) 



If the location xq is not known, the final state of the par- 
ticle will be described by a density operator proportional 
to 



p(x,x') 



f-aaPfao, f)f(x - X )*f(x' - Xq) dx 

S-^V^o, f) dx 



(2) 



This is properly normalized if f{x) is. We are implicitly 
assuming that the shape / of the localized wavepacket 
does not depend on the location xq of the localizing event, 
although the formalism could clearly be extended to deal 
with spatial inhomogeneities. 

At the risk of belaboring the obvious, we emphasize 
that the state described by J2J may still be "delocalized" 
in the sense that the probability p(x, x) to find the par- 
ticle at point x may still be nonzero over a large region 
of space (of the order of the original region covered by 
ij){x)). What has been mainly lost, in going from ip( x ) 
to the density operator (J5J, is the coherence between re- 
gions of space separated by a distance greater than I. As 
these regions can no longer interfere, it is legitimate to 
interpret @ as describing a situation in which the co- 
herent wavefunction of the particle is localized to a much 
smaller region of space of width /, only the precise loca- 
tion of the center of this wavepacket is not known with 
certainty. 

A simple example for which analytical calculations 
are possible is obtained if ip and / are both Gaussian 
wavepackets, ip(x) = (27r)~ 1//4 er 1 ' 2 exp(— x 2 /4(io) and 
f(x) = (27r) _1 / 4 /~ 1 / 2 exp(— x 2 /4l 2 ), where the original 
(Ax)o = cq may be much greater than I. We obtain, 
from Eq. QJ, 



p(xo,f) = 



2a l 
a[TT 2 



-x 2 /2(l 2 +a 2 ) = 



2<T l 



(3) 



defining the new a = (a 2 + I 2 ) 1 / 2 ; and, from Eq. J2J), 

1 1 tJ2 I „' 2 \ I At J* I j2\ 



p{x,x') = 



-(x 2 +x' 2 )/4(a 2 +l 2 ) 



/2tF y/a 2 + 1 2 

-(x-x') 2 a 2 /8l 2 (a 2 +l 2 ) 



x e 
e~ 



1 -(x 2 +x' 2 )/4* 2 -(x-x') 2 /8l 2 



r{x)^{x')e- {x - x '^' sl2 



(4) 



where the approximations apply in the case cfq 3> /. 
Equation J3J shows that, under these circumstances, the 
region where the particle may be found, as given by 
p(x, x), is substantially the same as before; the difference 
is in the spatial coherence represented by the off-diagonal 



terms (those with x ^ x r ), which is seen to decay with 
\x — x'\ over a region of the order of magnitude of the 
localizing length I. 

It is now easy to see that, in the limit in which 
I <C (Ai)o (where (Ax)o is the spread of the initial 
wavefunction) one should generically expect p(x, x') to 
take a form similar to The reason is that, in Eq. QJ, 
/(x—xq) will act as a S function, yielding a,p(xo, /) which 
is proportional to | -0 (a; ) | 2 (cf. Eq. Q); then, in Eq. J5J, 
p(xo,f) could be evaluated at xq — x (or xq — x 1 ) and 
factored out of the integral, which then becomes the au- 
tocorrelation of /, evaluated at x — x' . Thus, we expect, 
in general, 



p(x,x') 



\iP(x)\ 2 g(x-x') 
ip* (x)ip(x')g(x — x') 



(5) 



where g(x) is proportional to the autocorrelation of /. 
The second approximate equality is because we expect 
ip*(x)ip(x') ~ \4>{x)\ 2 if \x - x'\ < (Ax)o, which will be 
the case if \x — x'\ is restricted to a region of the order of 
I. 



B. Localization by momentum kicks 

Consider a particle with wavefunction ij){x) that re- 
ceives a momentum "kick" of magnitude fik at some point 
in time. The new wavefunction will be 



i>' k {x) = e ikx ^(x) 



(6) 



since the operator e tkx is the generator of a momentum 
translation. Assume that many different values of k are 
possible, distributed with a probability p(k). Then the 
density operator describing the state of the particle will 



be 



p(k)i// k (x)*i[>' k (x')dk 
P (k)e- tk{x - x } dk ^{x)*ip{x') (7) 



p{x,x') 



which is clearly of the form JSJ. If we search for the f(x) 
that p(k) is related to, we note that the Fourier transform 
of a convolution product is the product of the Fourier 
transforms, so p(k) ~ |/(fc)| 2 (where / is the Fourier 
transform of /). As p(k) is, by construction, positive, its 
square root is well-defined, and so one can simply take 



/(*) = 



1 



^ — ik 



x ^/pjk)dk 



(8) 



which is already properly normalized. For non- 
pathological p(k) , s, the function f(x) should have a width 
I ~ 1/Afc, so we can conclude that random momentum 
kicks lead to a "localization" of the coherence, in the 
sense of the previous section, within regions whose size 
I is inversely proportional to Afc, the range of momenta 
involved. 
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C. Representation of the localization process in a 
second-quantized formalism 

In order to apply these ideas to many-particle systems, 
especially bosonic systems, it is convenient to introduce a 
second-quantized formalism in which the possible states 
of each individual particle are represented by localized 
wavepackets, whose width is already of the order of mag- 
nitude of the width I associated with the localization 
processes. This "Wannier representation" approach may 
arise naturally in some important cases; in particular, for 
cold atoms in optical lattices. In a single-band descrip- 
tion, in which only one type of wavepacket needs to be 
considered at each site, the general pure state of a single 
particle can be written as 

IV) =^C i |0,...,l,...,0) (9) 

i 

with amplitudes Cj, where the "1" in the ket is at the i- 
th location and all the other entries are zeros. In this 
description, where the "lattice" is matched to the lo- 
calization process, it is easy to see what the action of 
a localization event is: a single basis state of the form 
|0, . . . , 1, . . . , 0) is left unchanged (the particle can only 
be localized at the i-th site in this case, and it is already 
localized there) , but the coherence of a superposition like 
(0 is totally destroyed, so p becomes 

p = ^|a| 2 |0,...,l,...,0)(0,...,l,...,0| (10) 

i 

We can now ask what happens when one has several par- 
ticles, of the bosonic type (so that more than one of them 
can be at the same site and with the same wavefunction) . 
We aim at describing the localization process using some 
appropriate "jump operator," with a view to eventually 
describing the dynamics of a system exposed to localiz- 
ing events by a master equation. The observation made 
above generalizes to the idea that pure Fock states of the 
form |m, n%, ■ ■ ■ , tin) must be left invariant by a localiza- 
tion event at the i-th location: the event can only happen 
if rii 7^ 0, but other than that the state does not change, 
since it already involves atoms whose precise locations 
are known. This means that the jump operator L, asso- 
ciated with a localizing event at site i must be a function 
of hi, the corresponding number occupation operator. It 
is also natural to make it proportional to hi, since this 
expresses the fact that a localization event is more likely 
to occur at a given site, the more particles one actually 
already has there. Note that this choice ensures that the 
single-particle complete loss of coherence, illustrated by 
the passage from © to (|10fl . still takes place: in that 
case, each jump operator hi can have only the eigenvalue 
or 1, so each one of them annihilates all the states ex- 
cept for one, and when the outcomes of all the possible 
events are added with equal a priori weights one recovers 
(|10|l . The novelty now is that a superposition such as, 



for instance, 

M = ^(|2,3> + |3,2» (11) 

is not rendered completely incoherent: a localizing event 
at site 1 transforms it into (2|2, 3) + 3|3, 2))/ vT3, and an 
event at site 2 transforms it into (3|2,3) + 2|3, 2))/vl3- 
Both events are just as likely, so the final density operator 
is 

p=i|2,3)(2,3| + i|3,2)<3,2| + ^|2,3)(3,2| + ^|3,2)<2,3| 

(12) 

This is no longer a pure state (since p 2 ^ p), so some co- 
herence has been lost, but overall, not much. Intuitively, 
this is because the localizing processes cannot very well 
distinguish between the state |2,3) and the state |3, 2). 
On the other hand, the state l|12[l is not, as will be- 
come clear in the following subsection, a terminal state of 
the evolution of the system: a sufficiently large number 
of successive localizing events would eventually destroy 
completely the coherence of all superpositions such as 
(|lip. Physically, this is consistent with the fact that the 
coherence in (jllf) is inseparable from the derealization 
of one of the five particles involved: we can say that, in 
the state lfTT|) . two particles are certainly at site 1, and 
two others at site 2, but there is a fifth particle that is, 
in some sense, in both locations at once. Localizing pro- 
cesses will eventually force it to settle in one site or the 
other, and at that point the coherence of the superposi- 
tion will be completely lost. From this perspective, the 
reason why it may take many localizing events to achieve 
this result is that there is a large probability that the par- 
ticles affected by the individual localizing events should 
be the ones that are localized already, as opposed to the 
elusive fifth one. 



D. Master equation treatment 

Having identified the jump operators for the localizing 
processes in the previous subsection, we can write down 
a master equation for the density operator of the second- 
quantized bosonic system on a lattice, in the standard 
Lindblad form 12], 

P = ~ r E ( L I L *P + P L l L * ~ 2L *P L l) 

i 

= -r*^2(h 2 p + ph 2 -2hiphi) (13) 

i 

where r is some appropriate rate, proportional to the 
probability per unit time that a localizing event may oc- 
cur anywhere. 

In subsequent sections (IV and V) it will be shown how 
equations of the form 113|) may, in fact, be derived, more 
or less from first principles, for some physical situations. 
For the moment, however, the focus will be on some of 
the formal consequences of (|13H . 
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To begin with, it is obvious that 1|13|) preserves any 
Fock state, as well as any incoherent superposition of 
Fock states. On the other hand, its effects on superposi- 
tions are easy to ascertain, because the hi operators do 
not mix different Fock states. Thus, for a general super- 
position given by a density operator of the form 

P(t) = ^2 Pn 1 ,...,n N ,n> 1 ,...,n' N (t)\ni, ■ ■ .,n N ){n[, . . . , n' N | 

(14) 

the equations for the matrix elements Pm,...,nx,n{,...,n' W 
all decouple, and we find 




(*) (15) 

Since (nf — n'^ — 2riin' i ) = (m — n£) 2 > 0, all the ampli- 
tudes Pn ll ... ! n Jv .n' 1 ,...,n' iv (t) eventually decay to zero except 
those for which, for every i, rii = n[, that is to say, the 
diagonal elements. Hence, as pointed out at the end of 
the previous subsection, all coherences between different 
Fock states are eventually destroyed, although some may 
decay much faster than others, depending on how large 
the differences (m — rt-) 2 are. 

E. Generalizations 

The most important generalization of the formalism is 
to deal with situations where the localized wavepacket 
may not be a perfect fit to the "natural" lattice for the 
problem. In the single-particle case, denoting by \i) the 
state |0, . . . , 1, . . . , 0) with a 1 in the i-th position, it is 
easy to see that the continuous treatment in Section II. A 
amounts, in the discrete case, to reducing a state of the 
form © to something of the form 

m=T,f(i'-m as) 

i' 

where /(£) is a discrete function of the (negative or pos- 
itive) integer argument i, centered at i = 0. The formal- 
ism developed in Sections II. C and II. D corresponds to 
the assumption that the particle's wavepacket collapses 
to a single lattice site, which is expressed by the choice 
f(i) = Sio (5ij is the Kronecker delta symbol). If the 
wavepacket spreads over a few lattice sites instead, one 
can get a result approximately equivalent to (|16|l by act- 
ing on the (single-particle) wavefunction with the opera- 
tor 

£i = £>'/(*'-*). (17) 

i' 

This produces f(i' - i)C v \i') ~ C, /(*' - 
provided the original state's coefficients C; do not change 
by much over the region where f(i' — i) is nonzero (cf. the 



remarks at the end of Section II. A and in Section II. B). 
The overall factor Ci is later removed by normalization. 

It seems, therefore, natural to generalize the concept of 
"localization-induced decoherence" to the many-particle 
case by adopting l|17[) for the Lindblad jump operators 
on the first line of Eq. Ijl^jl , and use the resulting master- 
equation to describe the dynamics when the localization 
processes cover a region somewhat larger than a single 
lattice site. 



III. DAMPING OF CENTER-OF-MASS 
MOTION BY LOCALIZING EVENTS. 

Consider a set of N identical particles that may occupy 
any of a total of M sites in a one-dimensional chain. 
The position of the center of mass of this system can be 
described by the operator 

M 

xcM = j-fJ2jhj (18) 
j=i 

if the lattice spacing is a. As is known from standard 
discretization schemes, such as the Bose-Hubbard model 
[l3L Il4| , an appropriate kinetic energy operator for this 
system is 

M 

K = -J^(a t j a j+1 + aj +1 & J y ) (19) 
i=i 

For convenience we shall assume periodic boundary con- 
ditions, so that M + 1 = 1 (i.e., all the additions of 
indices will be understood to be performed modulo M). 
The Heisenberg equations of motion, with the Hamilto- 
nian H = K, then yield the center of mass velocity 

dx C M «„, . , .Jay-^ fi t \ 

VCM = ~1T = h [ ' XCm] =l Nh^ [ a J+^ ~ a i a i+ l ) 

(20) 

If one applies [K,vcm] again, to calculate the rate of 
change of the center of mass velocity, one obtains a sum 
of terms that cancel identically, except for edge effects 
that are taken care of by the assumed periodic boundary 
conditions. This indicates that the Hamiltonian H = K 
is appropriate for a system of free (and non-interacting) 
particles. 

Imagine now that the system is somehow subjected 
to "localizing" processes, such as those considered in 
the previous section, that result in a master equation 
for the many-particle density operator p of the form 
given by l|13|) . Then we can calculate their effect on the 
time-derivative of the expectation value of any opera- 
tor O by the formula (d/dt)(0) = Tr(Op). The result 
for d{xcM) /dt is 0, since all the hi operators commute 
among themselves. On the other hand, the result for 
d{vcM)/dt is 
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d M 

-tt(ucm) = ( Tr ('"CMnfp) + Tr(v C Mpnf) - 2Tr(vcMriiphi)) 

i=l 



irJa 

' Nh 

2irJa 



Nh 



Af 



-2r(t) C M) (21) 



This means that localizing processes of the form consid- 
ered here result in a damping of the center-of-mass mo- 
tion of the system, at a rate (i.e., with a damping con- 
stant) 2r. To some extent, this should not be a surpris- 
ing result, since we showed in Section II. B that random 
momentum kicks tend to localize the particle (always in 
the sense of "breaking the long-range coherence of the 
wavefunction" ) in a way quite similar to actual random 
projective position measurements, and it has been known 
at least since Einstein's days that a particle subjected to 
random momentum kicks experiences an overall damping 
of its motion (fluctuation-dissipation theorem). What we 
have here is a formal derivation, and unification, of these 
various results, in a form suitable for use in quantum 
many-body physics. A specific example of localization 
by random momentum kicks is provided in the next Sec- 
tion, by a consideration of atomic spontaneous emission. 

Damping of the center of mass motion also results from 
the more general localization processes discussed in Sec- 
tion II. E, that is, processes that tend to localize the par- 
ticle over a region larger than one lattice site. With Li 
as given by l|17|l . direct calculation shows that the con- 
tribution of the term Lfp + pL\ — ILipLi to the time 
derivative of (aj +1 aj) is 

a ]+i a i L i + L 1 a ]+i a 3 - 2Liat +1 ajLiJ 

= (f(j ~ i) - f(j + 1 - *)) 2 (ot+iOi) ( 22 ) 

To calculate the total time derivative, one should sum 
(|22[) (and its complex conjugate) over all the values of 
i, which correspond to processes localizing the particles 
around the respective sites. When this is done, the pref- 
actor in ill' I'll becomes independent of j: 

i 

= E(/H)-/(l-*)) 2 ^/ 2 (23) 

i 

so that now 

i{voM) = -pr(v C M) (24) 



In general, less localization also means less damping. If 
f(i) — <5 i0 , one has / 2 = 2, as in Eq. (|2~T|l . On the 
other hand, if, for example, one takes the three-point 
function /(0) = l/y/2, /(±1) = 1/2, Eq. yields 
f 2 = l/2 + 2(l/ v / 2- 1/2) 2 = 0.586. 



IV. LOCALIZATION THROUGH 
SPONTANEOUS EMISSION (LIGHT 
SCATTERING) 

Consider a system of atoms in an optical lattice (as 
usual, we will consider one dimension only). If the laser 
beams used to form the lattice have wavelength A, then 
the lattice spacing is a = A/2. In general, the lasers 
are far detuned from the atomic transition, but it may 
still happen from time to time that an atom may ab- 
sorb a photon from the lattice and reemit it in a random 
direction. From the absorption, the atom receives a mo- 
mentum hk = ±27rh/\, whereas from the spontaneous 
emission the momentum component along the lattice di- 
rection can be anything between —2nh/X and 2nh/X. 
Overall, therefore, the atom experiences a random mo- 
mentum kick between —AnU/X and AnU/X. Most often, 
the whole process is referred to as spontaneous emission, 
but more properly it should be called "light scattering" 
(see |Tflj for a discussion of why this distinction may be 
important; see also f° r detailed descriptions of 

the decoherence induced by light scattering by free atoms 
or molecules). 

In a periodic lattice, the atom's momentum is quan- 
tized, and restricted to values 2irph/Ma, where M is 
the number of lattice sites and p is an integer between 
-M/2 and M/2 (first Brillouin zone). With a = A/2, the 
range of momenta from — Airh/X to Anh/X corresponds to 
~2Tr?i/a to 2TTh/a, twice the range of the first Brillouin 
zone. We shall take momenta that differ by 2nTi/a to 
be equivalent (i.e., quasimomcnta). For a single particle, 
we would describe the effect of a random spontaneous 
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emission event on a state of the form © as 

|^) = ^Q|0,...,1 ) ...,0) 

j 

^|V); = ^Qe 2 -^|0,...,l,...,0) (25) 



where the argument in the exponent is i(px)/h, as in 
Eq. 10, only x = ja for a particle at site j in the lattice. 

In the many-particle, second-quantized formalism, the 
natural generalization of (|25(l is to adopt a "jump oper- 
ator" of the form 0] 



,2iripj/M 



n ■, 



(26) 



This is most naturally seen by working temporarily in 
the momentum representation, where the creation and 
annihilation operators of an atom with momentum k are 
defined by 



1 



Ck 



M ^ J 



(27) 



Then if an atom in the momentum state k emits a photon, 
its momentum changes to k + p, and the appropriate 
operator is simply c^, +J9 Cfc. The sum of this over all k is 



M/2 ^ M/2 

M 



a -a,. 



k=-M/2 



E« 



Li 



k=-M/2 
,2wipj/M . 



(28) 



We can, therefore, use H2fi|l to write a master equation of 
the form 



i^(iJV + ^V2Vit) (29) 



p = —r 



where r is the (single-atom) rate of spontaneous emis- 
sion events, the sum over p is from —M/2 to M/2 — 1 
(since the quasimomentum M/2 is the same as —M/2), 
and the factor 1/M gives the probability that, in each 
event, the momentum kick received by the atom may be 
precisely p. This uniform probability distribution for p 
has been adopted for simplicity, as it immediately leads, 
via carrying out the sum over p, to 



(30) 



which is to say, the same as the simple master equation 
(|13[) derived in Section II. D. It follows that, under these 
conditions, spontaneous emission would cause damping 
of the center of mass motion with damping constant 2r, 
as in Section III. 



If the probability distribution g(p) of the momentum 
kicks is not uniform, one would obtain a somewhat dif- 
ferent master equation: 



P = —r g(i - j) [hihjp + phifij - 2hiphj) (31) 

where g(i—j) = ^2 p g(p)e 2nip ( l ~fi/ M is a discrete Fourier 
transform of g(p) . Note that this will be of the same form 
as the generalization discussed at the end of Section III, 
and based on the jump operator i|17[l . provided one has 
a function / such that 



(32) 



Ignoring possible complications arising from the finite 
size of the lattice, we see that Eq. (I32f) may be satisfied by 
taking f(J) to equal the Fourier transform of W g{p) (cf. 
Section II. D and Eq. JSJl). Thus, in this more general case 
also, and provided that the range of probable momentum 
kicks is sufficiently large, one obtains the effective local- 
ization of the wavefunction and the attendant damping 
of the center of mass motion, as discussed in Section III. 

The results in this Section are consistent with the in- 
tuitive notion that a spontaneously emitted photon car- 
ries information on the position of the atom, and there- 
fore tends to partially collapse its wavefunction. On the 
other hand, a simple picture of the atoms in a 1-D Bosc- 
Einstein condensate (such as the "tubes" to be discussed 
in the next section) as "billiard balls" might suggest that 
a spontaneous emission event simply gives a momentum 
kick that is transmitted down the chain of identical atoms 
and eventually leads to no other consequence that the 
ejection (from the condensate, if not from the actual 
physical lattice) of the last atom in the chain. Here we 
see that, in fact, the partial localization, and attendant 
decoherence of the many-particle wavefunction, caused 
by the emission event also results in an overall damping 
of the center of mass velocity for a condensate in motion. 

In a far-detuned optical lattice the "potential depth" 
V is proportional to fi 2 /A, where A is the detuning and 
fi the Rabi frequency, whereas the light scattering rate 
r is proportional to (Q/A) 2 T ~ VT/A, where 1/r is 
the lifetime of the excited state 01 • Since r is typically 
very small in current experiments, the motional damp- 
ing predicted here may not have been observed yet, but 
it should be a relatively easy matter to verify it experi- 
mentally, in a low-damping, pseudo 1-D geometry (such 
as the so-called "pancakes") by reducing the detuning 
A. If this is done while increasing Q, so as to keep the 
lattice depth V constant, one should observe a motional 
damping 7 ~ 1/A, regardless of other factors (such as 
the geometric considerations that might determine g(p) 
above). Note that this particular damping mechanism 
is, in fact, independent of the lattice depth (at constant 
scattering rate). 
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V. LOCALIZATION CAUSED BY 
RANDOMLY-FLUCTUATING POTENTIALS 

A. Formalism 

At least under some circumstances, random, fluctuat- 
ing potentials acting on a system of bosons may lead to 
master equations of the forms discussed in previous sec- 
tions. Consider a Hamiltonian with a term 



(33) 



The functions Vi represent varying, random potentials 
acting at different sites on the lattice. They will be taken 
to have zero average and be characterized by spatial and 
temporal correlations that will be more precisely consid- 
ered in what follows. 

The equation of motion for the density operator, 



P=-j-[H,p] 



(34) 



can be formally integrated once, over a short time interval 
[t — At, t], and iterated, to yield 



p = -i[H(t),p(t-At)}-^ 



t-At 



[H(tUH(t'),p(t')]]dt' 
(35) 

We are interested in the particular contributions to this 
equation arising from the term l|33[l . To this end we in- 
troduce a modified density operator which involves an 
average over different realizations of the stochastic pro- 
cess Vi. The first term in l|35l) then vanishes, because we 
do not expect the value of the external potential at t to 
be correlated in any particular way to the value of p at 
any time t' < t. This leaves 

1 

P^-pE / . ^(W(0>(*W(t') + P(t')fiihj 
hip{t')hj — hjp(t')fiij dt' (36) 



h — Jt-At 



Note that the above does involve a sort of Markov as- 
sumption, since p(t') and Vi(t) could certainly become 
correlated in general: for instance, p(t') could depend on 
Vi(t") at an earlier time t", and Vi(t) and Vi(t") could 
be correlated. Perhaps a better way to justify the ap- 
proximation is to simply assume that p(t) does not 
change very much over the time scale over which Vi(t) 
becomes uncorrelated. Under those conditions, we get 
the approximate equation of motion 

p ^ -7I yZ{ViVj) {hihjp + phihj - 2hiphj) (37) 
n 

hi 

where now everything is evaluated at time t, t c is a char- 
acteristic correlation time for the fluctuating potential, 
and (ViVj) is an equal-time average, i.e., we assume the 



space-time correlation function of the fluctuations factor- 
izes into a space part and a time part, as in, for instance, 

mwm = <V i V J > e -l«-* , IA*. 

Equation l|37|l is already formally reducible to cases 
considered earlier, at least for a homogeneous system (see 
Eq. (|31|l . in the previous section, and the discussion fol- 
lowing it), but for later purposes it is useful to consider 
explicitly what happens if we introduce at this point the 
momentum components of the fluctuating potentials Vi'. 
let 



M—l 



v, 



,2-Kijk/M 



(38) 



For a homogeneous system, we want (V 2 ) to be inde- 
pendent of j, and (VjVf) to depend only on j — j' . 
This means we must assume (VkVk>) = (\Vk\ 2 )Sk.-k' (it 
is understood that, because of the lattice periodicity, 
V-k = Vja-k)- Then we have 



2\ p 2ni(j-j')k/M 



(39) 



Introducing 



9ie 



-2-nikl/M _ 



-y 



* 2-nikl/M 
til e 



we find 



2\ e 2m{j-j')k/M 



M 



1,1' k 

= ^2gi9i'h',i-j+f 
1,1' 

= Y^9l-3'9i-j 
1 

1 

and l|37[) becomes an equation of the form 

P = -r J2 ( L \ L iP + P L \ L i - ZLipL] 
1 

with r = —T c /h 2 M, and the Lindblad operators 



(40) 

ikU-j'-l+l')/M 



Li =£*•- 



lUj 



(41) 



(42) 



(43) 



just like the ones discussed in Section II. E (compare 
Eq. (|17fl ). only the discrete function f(i) is now complex. 
The damping analysis in Eqs. I|22|) ~ (|23[) carries through, 
with the result 

J = r^2\g-i -gi-i\ 2 



/r * 

k 



h 2 M 



^E(^| 2 >) 1/2e ~ 27 " fe ' /M i 1 - e 2mk/M ) 



2\ • 2 

x sin 



nk 
~M 



(44) 
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in terms of the momentum components, 14, of the fluc- 
tuating potential V4 By undoing the Fourier transfor- 
mation (|38|) ■ this can also be rewritten in terms of the 
original Vf. 

i 

B. Application: damping of oscillations in 1-D 
Bose-Einstein condensates 

The results in the previous subsection may be used to 
provide a heuristic explanation of the damping observed 
in some recent experiments |9j on atomic Bose-Einstein 
condensates confined to an essentially one-dimensional 
geometry (a narrow "tube"), with a relatively weak pe- 
riodic lattice in the direction of the motion. It is known 
that the mean-field (Gross-Pitaevsky) description of a 
BEC in a 1-D lattice exhibits a dynamical instability 
when the wavenumber associated with the supercurrent 
flow becomes larger than 7r/2a (where a is the lattice 
spacing) Q. Although in the experiments of the 
average velocity of the condensate was kept well below 
the instability point, the large noncondensate fraction in 
these tightly-confined systems includes a non-negligible 
fraction of atoms with momenta larger than hir/2a. It 
has been argued by several groups l2fj| that this frac- 
tion of atoms is ultimately responsible for the damping 
observed in the experiments |2l[. 

A heuristic explanation might proceed as follows. The 
displacement of the lattice in the experiments turns the 
ground-state "vacuum" fluctuations into "real" excita- 
tions, without, however, substantially changing the frac- 
tion of atoms with momenta above the critical value 
(since the displacement is quite small). The chaotic mo- 
tion of this small fraction of atoms is seen by the rest 
as producing a random "external" potential, through 
the atom-atom interaction term (the term (U/2) J^. nf 
in the Bose-Hubbard Hamiltonian) . More specifically, a 
density fluctuation Srii at position i is seen as a fluc- 
tuating potential Vl = USrii by the bulk of the atoms 
at that site. We can, therefore, immediately apply the 
results of the previous subsection to conclude that the 
center-of-mass motion of the condensate will be damped 
as d{vcu)/dt = -2j'{vcm), with 

where r c is a characteristic correlation time for the fluc- 
tuations, and 5rik is given by the spatial Fourier trans- 
form of the fluctuations. (The factor of 2 in the damping 
7' = 7/2 has been introduced to facilitate comparison 
with the experimental results as well as with previous 
theoretical estimates.) As indicated at the end of the 
previous subsection, this can also be rewritten as 



The formula l|47(l is very similar to one introduced in 
lU, but appears to differ from it by a factor of 4, if one 
should simply identify the Srii above with the noncon- 
densate site densities hi considered in Q . This is mostly 
because in Q the operator nature of the fluctuation field 
was preserved for more steps in the calculation, until at 
the end certain products of bosonic operators were factor- 
ized using a standard ansatz that naturally brings about 
a factor of 2 (e.g., (ajaja.,- Oj + i) ~ 2(a]aj) (a]a J+ i)). This 
suggests that in trying to use (|46[) or H47[) for quantitative 
predictions, one may have to define the magnitude of the 
"external field" Srii in a somewhat ad hoc way. 

The range of momenta, Hit /2a < \p\ < hir/a, in the 
"mean- field unstable" region corresponds to M/4 < k < 
3M/4 in (|46|l . and this is centered on the maximum of 
the trigonometric function. This suggests that we could 
obtain an order-of-magnitude estimate of the expected 
damping by making the replacement 




where rihigh is the fraction of the N atoms in the tube 
with momenta higher than hi: /2a. This is a quantity that 
was calculated, for the experimental parameters, in 0. 
We may further assume that r c ~ ti/J, since J/h is the 
"hopping rate" between neighboring sites, and one may 
expect this to be the typical rate at which local density 
fluctuations would decay. This estimate, 

is compared in Figure 1 to the experimental data, with 
remarkably good agreement over a range of values of the 
lattice depth V where 7' changes by more than two or- 
ders of magnitude. (The formula used to fit the high- 
momentum data in Fig. 1 of 8] is rihigh = 0.01 + 0. 018 V+ 
0.0019y 2 ; see HQ for details on how to relate U and 
J to the lattice depth V . The total number of atoms 
N = 80 and the number of lattice sites M = 60.) 

The agreement shown in Fig. 1 suggests that this 
heuristic approach does capture the basic physics, at 
least in a qualitative way, for lattices that are not 
too deep. (For the deeper lattices, where the damp- 
ing seen in Fig. 1 appears to saturate, a very different 
pseudo-fermionization model developed in ]22\ is found 
to yield excellent quantitative agreement with the exper- 
iments.) A much more careful treatment of the present 
model, with a better factorization ansatz, and t c and 
{(Srii — 5ni + i) 2 ) calculated numerically from mean-field 
theory, can be found in |8j. 

It should be noted, however, that there are obvious 
limitations to the idea of treating the fluctuating field 
due to the atoms themselves as an equivalent "external 
field." As will be shown in the next section, the localiza- 
tion process described by a master equation such as i|13|) 
is equivalent to a sort of heating, that would eventually 
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1 2 3 4 5 6 7 



V/E R 

FIG. 1: Comparison of the experimentally observed damping 
(dots) and the estimate l|49fl (line), as a function of lattice 
depth V in units of the recoil energy Er. The vertical axis is 
in units of rad/s. 

destroy completely the coherence of the wavefunction on 
all length scales; in BEC terms, it would lead to total 
depletion of the condensate. In reality, however, the sys- 
tem of atoms in the experiment @ is, to a good approx- 
imation, closed and Hamiltonian, and the only external 
energy input it receives comes from the initial displace- 
ment of the trap. One would expect, therefore, that the 
modeling by an external fluctuating field would become 
inadequate once that small initial energy became "ther- 
malized" (if not earlier). 

VI. EFFECTS ON THE STATIC DENSITY 
DISTRIBUTION 

Although most of this paper has been concerned with 
the effects of a master equation of the form <|13f) on the 
system's dynamics, the processes leading to (fC^) could 
also have observable effects on the atomic density dis- 
tribution for a system at rest in a trap (inhomogeneous 
potential). This brief section takes a look at these effects. 

The Hamiltonian for a system in a trap potential given 
by the site function V(i) could be written in the form 

H = £ V{j)n 3 -JJ2 (ajfii+i + a] +l a 3 ) (50) 

3 3 

Now we are not assuming a finite number of sites nor 
periodic boundary conditions: the potential V(i) will de- 
termine which sites are, in fact, occupied. Direct in- 
teractions between the particles have been neglected for 
simplicity; a full treatment including interactions would 
be quite complicated and beyond the scope of the present 
paper. 

The full equation of motion for the density operator p 
will be 

^ = --^[H,p] - ( fl fp + P^i ~ 2 n t pn t ) (51) 

i 

This leads to a linear system of equations for the expec- 
tation values {a\aj) that make up the so-called "single- 



0.35 




FIG. 2: Evolution of the density distribution for a particle in 
a parabolic potential V(j) = Qj 2 with Q,/J = 0.1. Solid line: 
ground state density distribution, r — 0. Dashed line: density 
distribution after t — 4h/J, assuming r — 0.5J/R. Dotted 
line: a Gaussian fit to the central portion of the dashed line. 

particle density matrix" : 

j t {a\a 3 ) = -2r(l - 6 ij )(a\a j ) + i(V(i) - V(j)){&\&j) 

+ iJ (<4fi%) + (al-iOj) - (ojaj+i) - (ojfflj-i) 

(52) 

In the absence of the damping term r, the system would 
most likely be found in the ground state of the Hamilto- 
nian H , where all the expectation values (a\a 3 ) are pro- 
portional to the products C*C 3 of a single-particle wave- 
function of the form JjJ}. With r ^ 0, however, the sys- 
tem ()52|l becomes, formally, dissipative, and must evolve 
towards a final steady state in which all the (ffljflj) with 
i =/= j are zero, and all the (n 3 ) are equal (since the right- 
hand side of Eq. i|52|) involves the differences (hj±i)— (fij) 
when i = j± 1). That is to say, the ultimate steady state 
of (|52[1 for nonzero r is a "fully-depleted condensate," 
spatially completely homogeneous, even in the presence 
of a confining, inhomogeneous potential 

The explanation for this last result must be sought in 
the fact that the processes responsible for fl^jl represent 
a sort of heating, and if they are allowed to continue in- 
definitely they will eventually provide the system with 
enough energy to render the confinement potential V ir- 
relevant. In practice, one expects that the description 
(|13fl will cease being even approximately valid long be- 
fore a totally homogeneous state is reached; but for short 
times (assuming that one could "turn on" the r term in 
Eqs. (|52|l at will, by, for instance, reducing the detun- 
ing and bringing up the spontaneous emission rate, as 
discussed in Section IV), the distinctive effect of the pro- 
cesses (|13|l on an inhomogeneous system at rest will be a 
flattening of the spatial distribution. 

This effect is illustrated in Fig. 2 for a system of atoms 
in a parabolic trap. The potential is chosen so that the 
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initial ground state (with r = 0) is given by the solid line, 
which is approximately Gaussian (as one would expect 
for the ground state of a harmonic oscillator). Then the 
r term is turned on, and after a time equal to 8/r the 
new density distribution is given by the dashed line. To 
see that the effect is truly a flattening, and not just a 
broadening, a Gaussian with the same curvature at the 
center has also been plotted, for reference (dotted curve): 
it clearly is much broader than the actual solution. 

There are a number of ways in which this effect might 
be useful. For instance, it might be possible, by watching 
the shape of the distribution change, to estimate the rate 
at which "localizing processes" of the type considered 
here are taking place in a given system. Another intrigu- 
ing possibility might be to deliberately, and temporarily, 
"turn up" the value of r before driving a trapped BEC 
through the Mott insulator transition, and see if this re- 
sults (necessarily after further cooling) in a more uniform 
filling of the lattice in the insulator phase — a sort of an- 
nealing approach. More detailed calculations, including 
interaction terms, would be necessary to really assess this 
possibility. 

VII. CONCLUSIONS 

In this paper, a formalism has been developed to de- 
scribe the loss of long-range coherence of a many-particle 
wavefunction due to random "localizing" events, includ- 
ing the physically relevant case of random momentum 
kicks. The resulting master equation formalism is espe- 



cially appropriate to the study of cold atoms in optical 
lattices. It has been shown that these random localiz- 
ing events lead to a damping of the center of mass mo- 
tion of the system, and several types of processes have 
been identified as producing these effects: in particular, 
spontaneous emission (or, more precisely, light scatter- 
ing), and random external potentials with the appropri- 
ate correlations. The latter have been used to provide a 
heuristic model of the recently-observed strong damping 
of the motion of one-dimensional BECs in a relatively 
weak lattice. Finally, the signature that these processes 
might have on the density distribution of a condensate 
at rest has also been considered. 

The relatively simple calculations presented here could 
clearly be extended in several ways, in order to deal with 
more realistic systems: for instance, using the detailed 
momentum distribution of scattered photons in spon- 
taneous emission, specific examples of random external 
potentials, the effects of interactions. . . If the formalism 
introduced here is found to be sufficiently promising for 
dealing with real- world problems, these issues will almost 
certainly be explored in further work. 
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